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Abstract 

Absorption of light by a nanoparticle in the presence of resonant atom 
and fluorescence of the latter are theoretically investigated. It is shown, 
that absorption of light by a nanoparticle can be increased by several or- 
ders because of presence of atom. It is established, that optical bistability 
in such system is possible. 

PACS numbers: 42.50.Ct, 12.20.-m, 42.60. Da, 42.50.Lc, 42.50.Pq, 42.50.Nn 

1 Introduction 

The cross-section of light adsorption by an isolated spherical nanoparticle imbed- 
ded in a host medium and which radius a is essentially smaller then light wave- 
length A in the medium (a/A 1) is given by the classical formula 1 

where S = ira 2 , e = e p /eh = e' + ie" is the relative complex dielectric function 
of the nanoparticle, Sh and e p are dielectric functions of the nonabsorbing host 
medium and nanoparticle respectively. 

As a rule, a p is smaller then geometrical cross-section of the nanoparticle S. 
On the other hand, it is well known that the cross-section of resonant atom-light 
interaction is considerably larger: 

Z7T 7 

where 70^ and 7 are the radiation and the total width of the resonant transition 
of an atom imbedded in a host medium. Note that 70?! is expressed in terms 
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of free-space spontaneous emission rate 270 as 70/1 

e.g., mw 

As a rule, a p <C a a (see Fig.[Q. 



(e/ l ) 1/2 |(e/ l + 2)/3|7o (see, 




Figure 1: Atomic lens. 



The aim of the paper is clarification of the probability of cascade energy 
transfer from light to an atom and then to a nanoparticle 0], and investigation 
of atomic fluorescence in this conditions. 

2 Oscillating classical dipole 

Let us consider an auxiliary problem connected with atomic excitation transfer 
to a nanoparticle, the particle absorption of the electric field energy of clas- 
sical dipole d = d{sin(ip),0,cos(ip)} which oscillates with frequency uj and is 
at a distance R from the nanoparticle center. The power Q c absorbed by the 
nanoparticle can be represented in the form 



where <5E(R) is the 'image field' of the dipole, and overline denotes time aver- 
aging over the time that is considerably greater than period of light wave. Since 
<5E depends linearly on d, the absorbed power Q c can be rewritten as 



Q. 



d(R) <5E(R) , 



(3) 
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where G Q( g(R, R; ui) is the field susceptibility (or tensor- valued Green function): 

SE a (r, u)=J2 G *p( r > r '; ")M*', w) • (5) 





2.1 Field susceptibility 

In the near zone (R <C A) the 'image field' of the dipole (and, consequently, G) 
can be found by solving electrostatic Laplace equation. One can start with the 
scalar potential of a single charge e located on the axis Z at the distance R from 
the center of the particle 0): 



p e (r,R) = --£ 
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where P n (cos(#)) is the Legendre polynom and 6 is an elevation angle of the 
vector r. 

The potential of the point-like dipole d located on the axis Z at the distance 
R from the center of the particle is the sum of potentials JSJ) caused to nearly 
situated charges e and — e 
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where is an azimuth angle of the vector r. 

So, the 'image field' and G are found from JJJ and represented as a series 
(see, for example, This series diverges when R — > a, so that the higher 

terms start play the major part in it. Fortunately, one possible to rewrite it in 
a reasonable way, so that the field susceptibility tensor can be expressed in the 
form 
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where p = e/(e + 1), £ = pp , p = (p - p x ), p = i?/a, v = l/(e + 1), 
e = Ep/eh, £ P and are complex dielectric function of the nanoparticle and 

a, b 
c 

Let us concider some limit cases. 



host surroundings accordingly, F 



s is hypergeometric function. 
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|e| ^> 1 (i.e., ideal conductor) 



G ZZ (R) 



1 
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(10) 
(11) 



As is well known, in this case the 'image field' that corresponds to 110(1 and 
((11(1 can be represented by the sum of the fields of the charge go = ad z /R 2 
and — g , and of the dipole d = (a/R) 3 d{— sin(^), 0, cos(V>)} field. The 
charge -go is placed at the particle center; the charge go and dipole do are 
located along Z axis at the distance ag from dipole d(another words, at 
distance a 2 /R from the particle center toward the dipole d). 

Similarly to this limit case, the first three terms in square brackets of the equa- 
tions JSJ) and @ can be interpreted as the sum of the fields of the charges 



A* (e ^ l)/( £ + l)?o, -q, g e 



-q/e, —q e , and of the dipole di = /i J (e 



l)/(e+ l)do. Charge q and dipole di are located along Z axis at the distance 
R[ = ag from dipole d, whereas the charge — q is disposed at the distance 

1/2 

Ri = ap [g/(p/fi — p^ 1 )] from dipole d along the same Z axis. Charges g e 
and — g e are located at the distance R\ and R[ respectively from dipole d along 
X axis (see Fig- El) - These terms of the equations (JSJ and 10 dominate while as 




Figure 2: Geometry of the problem. 

R/a > 1, and e is far away from the region of the surface multipolc resonances 
that take place when Re(n + 1 + v) = 0. 
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Small distances, R/a — 1 <C 1 (go — 0) 



fV, , (R) ^ i G ZZ (R) ~ ^4 4' ( 12 ) 



2 K ' e+1 e h R 3 g 
It is exactly the case of the planar interface 
• Large distances, fi/a>l (i.e., g ~ R/a) 



G xx (R) ~ i G zz (R) ~ - — - -^-r . (13) 

The 'image field' corresponding to (|13fl is given by 

5E(R)^i^-^[3(nd)n + d] , (14) 

where n is the unit vector in the direction from the center of the particle 
to dipole d. 

Physical interpretation of the formula (|14H . At large distance the electric 
field Ed of the dipole d is homogeneous in the vicinity of the particle. The field 
of the polarization (or scattered field) <5E of the particle in such homogeneous 
field Ed at the location of the dipole d is given by (see, for example, £Q) 

<5E(R) = i^i iL [3(nE d )n - E d ] . (15) 
In turn, the field Ed in the quasistatic approximation is given, as is well known, 

by[U 

E d = -^3 [3(nd)n - d] . (16) 
Substituting this expression in l(T5|) we get {T3J|. 



3 Quantum consideration 

Energy transfer to the nanoparticle from the real atom exited by the light is 
given by (instead of (@J) 



Q c = 2lj Im 



^G Q/3 (R;c)(: o^d+ :) 

a, /3 



(17) 



where symbols ': : ' and '(•••)' denote the normal ordering operator and the 
quantum averaging respectively. For two-level atom d + is 



i+ - 



= d nm exp (itt t) a. 



(18) 
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where fig = u> — u> mn , d mn and u> mn are matrix element of the dipole moment 
of the m — n atomic transition and the resonance frequency of this transition, 
a± are the raising and lowering Pauli's operators. 
Substituting l(T%|) in l(T7|l results in 



Q c = 2lo Im 



^G a0 (R; W )d 



a, (3 



(19) 



where p m is population of the upper atomic level m, and j c is the addition 
nonradiative broadening of the atomic transition due to energy transfer from 
the atom to the particle (see below). 

4 Density matrix 

Density matrix for two-level atom in the vicinity of the nanoparticle obey the 
follow system of the equations [7] 



^ = -( 7 -ifi)p-in B A/2, 

^ = -2 7 (l + A)-2Im[iW] . 
at 



(20) 
(21) 



where A = p m — p n = (03) and p = (a-) exp(iu/t) are the population difference 
and coherence of combining levels, £Ir — |E • d mn \/h is Rabi frequency, E is an 
local electric field acting on the atom, Q — fl + v, v = Re(r), 7 = Im(r), 



a,/3 



(22) 



Here G ex is the exact field susceptibility. Therefore the imaginary part of Eq. 
1)22(1 as the total decay rate, 7 = 7r + j c , describes both radiative decay 7 r and 
nonradiative one j c - 

The field susceptibility, G, Eqs. is the quasistatic approximation 

of the exact one, G ex . Hence, it is responsible for nonradiative decay only. 
Nevertherless, the radiative part of the decay rate can be found in our case 
by use of quasistatic solution 10. Indeed, the solution represents multipole 
expansion of the scalar potential induced by dipole d. In the near zone the 
dipole part of the induced potential JJJ) is described by its term with n = 1 



p« (r, R) = ^ -|g [2 cos(V>) cos(0) - sin(^) sin(0) cos(0)] = %^ 



a ^ e ~ 1 a 



--{-sin(V0,0,2cos(V0}. 



(23) 
(24) 
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Dipolcs d and d p oscillate in phase due to the inequality a, R <C A. So, the 
emission probability and intensity of radiation are proportional to the total 
dipole squared |d + d p | 2 . Hence, the radiative part of the spontaneous decay 
rate of an atom placed next to a nanoparticle is given by 



Tr = Joh 



1 



1 a 6 



2 R 3 



1 



1 a 3 



2 R 3 



cos 2 (V>) 



(25) 



Expression (|25J) agrees with the radiative part of the spontaneous decay rate 
calculated without assumption of the quasistatic approximation in the limit 
a/A ^0|. 

Steady-state solution of the equations l|20|) . I|21|) is conveniently expressed 

i rvmm (26) 

where following dimcnsionlcss quantities are introduced: the total broadening 
of the transition 7 = "f r /"f h + Jc/^oh 



e-1 a 3 
~ e + 2 R 3 



e-1 a 3 

1 + 2 ^ 

e + 2 R 3 



cos 2 (V>) 



k- 3 Im [G ZZ (R) cos 2 (^) + G XX (R) sin 2 (?/>)] 



(27) 



(k is the wave vector of the light in the medium, tJj is the angle between R and 
d ~ E); local field gain factor £ = \E/E \ 



e + 



1 a 6 



2R 3 



[3(e ■ n)n - e ] 



(28) 



(Eq and eo are amplitude and unit polarization vector of the incident light wave); 
formfactor of the optical transition line g(fl)/ir = vr -1 / (1+U 2 ); = fi/(7 r +7 c ); 
and intensity of the incident light / = I/I s , where I s = fiw"foh/o' a is saturated 
intensity of free atom in the medium, a a is resonance cross-section (J2Jl. The 
second term in the braces of (I28|) multiplied by Eq is the scattered field of the 
nanoparticle in the near zone. 



5 Efficiency of the cascade energy transfer 

Efficiency of the cascade energy transfer 77 = Q C /Q P (where Q p — Ia p ) can be 
found by using expressions (|TT)|l and l(2*B|) 

1 a <7 a | £ +2| 2 \i c r 2 e 9 m ( 

V 24ir a S e"I l + ^^g(n)l' [ ' 

Obvious, efficiency of the cascade energy transfer can rich as much as several 
orders of magnitude because of multiplication of large values A/a, a a /S, \e + 
2\ 2 /e". 
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The limit cases 



• Large distances, R/a ^> 1 

In this case 7 — > 1, % ~ R~ 6 (see Eqs. (J27J and (H2J). Hence, Eq. is 
became 

77 ~ iT 6 . 

• Small distances, JJ/o - 1 « 1 

In this conditions 7 « 7 C — (_R/a - a/R)~ 3 (see Eqs. |2?J| and l(TT)l1, 
Assuming in addition i7~ 2 £ 2 g(f2)I <C 1, we get from Eq. 

77- (R/a-a/Rf. 

This decreasing of 77 is rather unexpected because the probability of the 
energy transfer from atom to nanoparticle is about 1 in this case. How- 
ever, atomic resonance cross-section is decreased when 7 is increased (see 
Eq.©). 




Figure 3: Efficiency of the cascade energy transfer as a function of the dimen- 
sionless distance R/a between resonance atom and centre of the nanoparticle. 
It is supposed that the light frequency is tuned in to the resonance atomic 
frequency at any distance. Calculations are made for a silver nanoparticle 
(e p = -15.37 + £0.231, A = 6328 nm) 

So, the efficiency is decreased both for large and small distances between 
atom and particle. Therefore, it achieves a maximum at an intermediate R > 
a. Fig. |3 shows efficiency of the cascade energy transfer 77 versus R/a in the 
assumption of the exact resonance at any R/a. On the other hand, figure Fig.Q] 
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Figure 4: The same as in Fig.|3| but for fixed light frequency 17 = 17/ 70 = —410. 

demonstrates very sharp resonance dependence of rj as a function of R/a when 
the frequency of the light wave is fixed. 

This sharply outlined resonance can be used to determine location of an 
atom regarding the surface with subnanometer precision. 

6 Fluorescence 

As it is known, the intensity of fluorescence If is proportional to "f r p m . Using 
Eqs. (|25|l - lt2fi|) . intensity of fluorescence can be written by 



(30) 



Fig.[S]shows the intensity of fluorescence, a.u., as a function of R/a. Conditions 
are the same as in Fig. 0] 



7 Heating of the particle 

Basic approaches and approximations 

• Steady-state approximation 

• Uniform temperature T p inside the particle 

• Laplace's equation Al), = for the temperature of the host medium 
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Figure 5: Intensity of fluorescence, a.u., as a function of R/a. Conditions are 
the same as in Fig. 0J 



• Energy balance equation 
Q c + Q P = x J VT h dS 

where x is thermal conductivity of surroundings 

• Linear temperature dependence of e p : 
4 = <o + a ( T P - To) 

• \e';/(e' p + e' h )\, \e h /e> p \^l 

Solution of the Laplace's equations is 

T h (r)=T + ^(T p -T ). (31) 
It gives a linear temperature dependence of heat removing from the particle 

Qt = Qc + Q p = ^xa{Tp - T ) • (32) 

Substituting (JQ and (|29|l in this equation results in following cubic equation 
with respect to the relative increase of the image part of the dielectric function 
of the nanoparticle z = 8e' p / e' pQ : 

(x-y-l)z 3 + [3(x -y)-2] + [3(x - y) - 1 + (x - 1)/] z+xf+x-y = , (33) 

where the following dimensionless quantities are introduced: / = i^~ 2 ^ 2 ^(fl) J 
is saturation factor, x = [N p +l)z/N, y = Az/N, z/N = [fru) mn 'y c / (Anax)] [a/e'' ], 
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Npo = Q P o/ ' (fu^mnlc) and N = Qj- / (Hu) mn 'y c ) are the number of photons ab- 
sorbed by the nanoparticle during the time 7J 1 directly from the light wave and 
the total one respectively. 

As well known, the qubic equation ll-i-it may have 3 solution at some param- 
eters. Therefore, z may exhibit bistable behaviour. In the Fig. 6 it is shown the 
regions of such bistability. 

Bistability Regions 



0.65 
0.6 

0.55 
0.5 



12 14 16 18 20 22 24 

f~I 

Figure 6: Bistability regions of the the relative increase of the image part of 
the dielectric function of the nanoparticle z = 5e p /e p0 , The range of values 
of parameter y = Az/N which depends on x = (N p + \)z/N and saturation 
factor / = ^f~ 2 ti 2 g(£l)I is different in each domain (N p q = Q p o/(huj mn j c ) and 
N = Qt I (fudmnjc) are the number of photons absorbed by the nanoparticle 
during the time 7^ directly from the light wave and the total one respectively). 
However, parameter y is always negative in all domains of bistability. 



8 Conclusions 

• Cascade energy transfer efficiency can rich as mach as several order of 
magnitude (10 3 - 10 5 ) 

• Efficiency is drastically decreased at both large and small distances be- 
tween atom and nanoparticle surface 

• For constant light frequency the efficiency as sharply as resonance depends 
from the distance between atom and nanoparticle 

• This sharp dependence can be used to determine the atom position near 
the surface 

• Bistability may take place when the population difference and the rela- 
tive growth of the image part of the particle dielectric function have the 
opposite signs. 
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